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The Raman response in the antiferromagnetic 2-leg S = 
1/2 Heisenberg ladder is calculated for various couplings by 
continuous unitary transformations. For leg couplings above 
80% of the rung couphng a characteristic 2-peak structure 
occurs with a point of zero intensity within the continuum. 
Experimental data for CaV205 and LayCai4_yCu2404i are 
analyzed and the coupling constants are determined. Evi- 
dence is found that the Heisenberg model is not sufficient to 
describe cuprate ladders. We argue that a cyclic exchange 
term is the appropriate extension. 

PACS numbers: 75.40. Gb, 75.50.Ee, 75.10.Jm 



Strongly correlated electron systems in low dimensions 
are of fundamental interest due to their fascinating prop- 
erties resulting from strong quantum fluctuations [0-p|. 
Important experimental insight is gained from spectro- 
scopic measurements of such systems. The spectral den- 
sities measured yield information on the kinetics and on 
the interaction of the elementary excitations as well as 
on the matrix elements involved. Thus quantitative the- 
oretical calculations of spectral densities are a major task 
in condensed matter physics. We use optimally chosen 
continuous unitary transformations (CUT) to map com- 
plex many-body problems to a tractable few-body prob- 
lems Q. This clear concept serves as a perfect basis to 
compute spectral densities of strongly correlated systems 
thus establishing a quantitative contact between theory 
and experiment Q. 

We will focus on optical investigations, in particu- 
lar on the Raman response, of antiferromagnetic 2-leg 
Heisenberg ladders realizing quasi one-dimensional (ID) 
strongly correlated systems. There are several experi- 
mental realizations of spin ladders like CaV2 05, SrCu2 03 
and LayCai4„yCu2404i rendering direct comparison be- 
tween theory and experiment possible ||7Hll| ■ 

Raman scattering measures excitations with zero 
change of spin and momentum. Starting at T = from 
the 5 = ground state the singlet excitations at zero mo- 
mentum are probed. The Raman response in spin ladders 
was recently calculated by first order perturbation theory 
for spin ladders |]12[ and by exact diagonalization ||l^ . In 
this work, we present detailed predictions obtained from 
CUTs using rung triplets as elementary excitations. Our 
results are not resolution limited because neither finite 
size effects occur nor artificial broadenings are necessary. 

The Hamiltonian for the 2-leg Heisenberg ladder reads 



where J|| > and Jj^ > are the leg and rung couplings; 
the subscript i denotes the rungs and 1, 2 the two legs. At 
T = the Raman response /(w) is given by the retarded 
resolvent 



I{uj) = -TT^^Im (0 \R\uj - H + iS)-^R\ O) 
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The observables i?'"'^^ (i?'°s) for magnetic light scatter- 
ing in rung-rung (leg-leg) polarization read in leading or- 
der |l|,|l| 
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The factors Aq^ and ^™"^ depend on the underlying 
microscopic electronic model. It is beyond the scope of 
the present work to compute them. Equally, we do not 
consider resonating Raman excitation processes. Results 
will be given in units of the factors squared. 

Technically, we employ a CUT to map the Hamilto- 
nian H to an effective Hamiltonian Hcs which conserves 
the number of rung-triplets, i.e. [iJoj^cff] = where 
Hq := H\j..^o IqmI- The ground state of iJoff is the 
rung-triplet vacuum. For the response function I{lu) the 
observable R is mapped by the same unitary transforma- 
tion as the Hamiltonian to an effective observable i?cff- 
We implemented the CUT perturbatively in x := J\\/J± 
and calculated Heft to high orders (1-triplet terms: 14'^, 
2-triplet terms: 13**^ order). The effective observable R^s 
is computed to order 10 in the 2-triplet sector and to or- 
der 7 in the 4-triplet sector. Generally, higher orders 
make higher accuracy possible. As a rule of thumb, a 
calculation in order n accounts for hopping or interac- 
tion processes extending over a distance of n rungs. 

The truncated series gives quantitative results up to 
X fa 0.6. Using standard extrapolation techniques like 
Pade approximants and optimized perturbation theory 
||l7| , |l8| , the effective operators R^s and H^ff can be cal- 
culated up to X « 1 with an uncertainty of about 5%. 
A qualitative description is obtained for x « 1.2. The 
Raman spectral density is calculated as continued frac- 
tion by tridiagonalization [[l9| . Because in ID asymp- 
totically free particles with quadratic dispersion display 
square root behavior at the band edges (van-Hove singu- 
larities) we use a square root terminator for the continued 
fraction. Thus neither finite size nor finite resolution af- 
fects our results. No divergences occur because a rung 
cannot be excited twice, i.e. triplets exclude one another. 



The relative motion of a pair of triplets corresponds to 
the dynamics of a particle on a half-infinite chain [|Tj . 

Sectors with odd number of triplets are inaccessible by 
Raman scattering due to the invariance of the two ob- 
servables R^g and -Rg^"^ with respect to reflection about 
the centerline of the ladder. Thus only excitations with 
even number of triplets matter. Therefore the leading 
contributions to the Raman response come from the 2- 
triplet and the 4-triplet sector. The total spectral weight 
(integrated over frequencies and momenta) of these con- 
tributions is depicted in Fig. n^ mn. The 4-triplet sector 
has a spectral weight I4 of less than 7% of the spectral 
weight I2 of the 2-triplet sector at a; « 1. Hence we focus 
on the 2-triplet contribution in the sequel. 
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FIG. 1. Spectral weights /„ as function of a; = Jy /J± ; n de- 
notes the number of triplets involved; the superscript denote 
the observable as defined in Eq. (H). 



In Fig. g the spectral density I{uj) from the 2-triplet 
sector is shown for various values of x = J\\/J±- The 
line shape is the same for i?''""g and R^°^ pfl] because the 
Hamiltonian is a weighted sum of the two observables 
H ^ i?™"s + xR^"^ (for Ao = l). Thus the excited state 
jjrungjQ^ equals — xi?'°^|0) except for a component pro- 
portional to the ground state |0) which does not matter 
at finite frequencies. This fact leads also to the intersec- 
tion at X = 1 visible in FigJ^. 

The spread of the lines in Fig. 2 on increasing x in- 
dicates clearly the increasing band width. For small x 
the Raman intensity shows a strong resonance near the 
lower band edge. This resonance is a consequence of the 
2-triplet attraction on neighboring sites pl^,&31 . Above a 
certain finite total momentum this leads to a 2-triplet 
bound state [2l]-p6lO| of which the resonance at the 
lower band edge is a precursor. Since for larger x the 
kinetic energy of the relative motion of the triplets in- 
creases the influence of the attraction decreases. There- 
fore, the resonance is rapidly broadened and shifted to 
the center of the continuum. In view of analyses of the 
spin gap ||l^ we note that it is not possible to detect the 
onset of the 2-triplet, non-resonant Raman continuum, 
i.e. twice the spin gap, reliably for x ^ 0.4. Furthermore, 
we found that the non-resonant line shapes do not de- 
pend very much on the precise form of the excited state 
R\Q). The qualitative features depend more on kinetics 
and on interaction. 




FIG. 2. Raman spectral density I{uj) for various coupling 
ratios x. For each curve the total weight is set to 0.1[Ao]. 
Long dashed lines depict lower and upper band edge of the 
2-triplet continuum. The short dashed line shows 2Lu{k = 0) 
for X > 0.5 where this is not the upper band edge {u-i(k) is the 
1-triplet dispersion, see inset in Fig. W). 



In Fig. 2 for x > 0.6, a second peak is visible near 
the upper boundary of the 2-triplet continuum, becom- 
ing more pronounced on increasing x. This feature is the 
combined effect of 1-triplet kinetics, 2-triplet interaction 
and matrix elements. First, the occurrence of a dip in 
the 1-triplet dispersion w(fc) at fc = (cf. inset in Fig. ^) 
leads to an additional van-Hove singularity situated at 
2u;(0) providing additional spectral weight. The impor- 
tance of this effect is illustrated by the short dashed line 
in Fig. 0. Second, the additional spectral weight is sepa- 
rated from the main peak by a double zero in the spectral 
density, see Fig. g. This double zero stems from a simple 
zero in the matrix elements implying that at a certain 
energy oj the state i?|0) is orthogonal to the excited state 
\uj) ! This intriguing phenomenon results from destruc- 
tive interference between several coupling contributions. 
We found that the destructive interference is triggered 
by the 2-triplet interaction since it vanishes when the 2- 
triplet interaction is switched off by hand |2^. Hence 
the orthogonality is induced by the interaction recalling 
in a broad sense an orthogonality catastrophe. Indeed, 
an arbitrarily small amount of the interaction suffices to 
induce at least a very narrow dip with the double zero at 
its bottom. We are led to the conclusion that the large 
density of states provided by the additional van-Hove sin- 
gularity renders the system particularly susceptible to the 
interaction- induced orthogonality. 

Analyzing experimental data, in Fig. |3| the Raman shift 
for CaV205 M is shown. This substance is a quasi- 
2D layered material where the S = 1/2 V^"'' ions form 
weakly coupled 2-leg ladders (trellis lattice). Susceptibil- 
ity measurements 128] predict weakly interacting rungs 
with X « 0.1. Theoretical analysis in first order yields 
X — 0.11 and J = 447cm^^ Il3|. For these values of x 



our results can be considered to be exact. It is instruc- 
tive to compare fits based on our results to the previ- 
ous analysis. Two fits and their parameters are given in 
Fig. |3[ Assuming the same experimental resolution as 
in Ref. ^ we find x = 0.09 deviating by 20% from the 
value obtained in first order. A change Ax of order x^ 
was to be expected since terms of order x'^ were neglected 
in the first order analysis |1^ . Hence the relative change 
Ax/x is of the order of x. We have to stress that the fits 
are extremely sensitive to the experimental resolution as- 
sumed. Assuming zero resolution (F = 0, see Fig. 0) the 
best value of x is 0.125, i.e. it is changed by about 40%. 
This remarkable sensitivity results from the very narrow 
resonance the height of which changes quickly as func- 
tion of the resolution. An effective value of F adding to 
the experimental resolution might be generated by resid- 
ual couplings not considered explicitly. It is unfortunate 
that the dominant phonon at 940cm~^ prevents the ob- 
servation of the upper part of the continuum so that the 
theoretically possible high precision analysis cannot be 
performed. Within the non-resonant Raman theory, our 
results exclude that the shoulder at 970cm~^ besides the 
dominant phonon is of magnetic origin. We conclude that 
the previous analyses Pq,pUl3jl2] and ours agree within 
the achievable accuracy. 



or extensions of the model are necessary. 
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3. Raman response for CaV205 from Ref. 



Solid 



mental resolution of F = 3cm . Long dashed line: fit with 
X = 0.125, Jx = 447cm~\F = 0. Thin constant line: offset 
of the fits to account for background. 



In Fig. ||the Raman lines [Q for La6Ca8Cu2404i are 
analysed. La6Ca8Cu2404i is a layered material contain- 
ing Cu02 ID spin chains and CU2O3 2-leg spin ladders 
p9[ . The inter-ladder coupling is weak and frustrated 
(trellis lattice) so that the ladders can be treated as 
isolated ladders. Because the atomic distance between 
neighboring copper sites in rung and leg direction is al- 
most the same one expects the spin ladders to be in the 
isotropic regime x ~ 1. This view is corroborated by 
the analysis of the 2-triplct bound state observed by in- 
frared (IR) absorption |0] leading to values of J± be- 
tween 1020 and 1100 cm . On the other hand, the spin 
gap values and neutron scattering data cannot be recon- 
ciled with the model in Eq. |] for J|| « Jj_, see e.g. Q 
and references therein. Either values of x beyond unity 
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FIG. 4. Raman response for La6Ca8Cu24 04i. Rough 
thin lines: experimental data ||lO| in (aa) (dashed) and 
(cc) (solid) polarization scaled to the same constant value 
between 2000 and 2500 cm~^ . Thick solid line: theory 
for J|| — J± — llOOcm^^ and resolution F — 3cm~^. 
Thin constant line: offset of the fit to account for back- 
ground. Inset: 1-triplet dispersions u}{k) for Jy = J± and 
a^cyc ~ 0,0.05,0.10,0.15 (solid, dashed, long dashed, dotted). 



Assuming the model in Eq. n and the observables in 
Eq. Is] the line shapes for the two polarizations should 
be identical independent of the value of x as explained 
above. The fact that this is not the case, see Fig. ^, indi- 
cates that an extension of the model is necessary though 
it cannot be excluded that a modification of the observ- 
ables (0) would also explain the deviations. 

The X = 1 result agrees quite well with the experi- 
mental data. The value J± — llOOcm^^ is consistent 
with the IR result |lj]. The main peak is situated at 
the right energy and has approximately the right width. 
The weight of the experimental high energy shoulder cor- 
responds to the weight in the theoretical second peak at 
higher energies. However, important deviations remain. 
(i) Experimentally, there is a high energy shoulder but no 
second peak, (ii) The main experimental peak is sharper. 
Inspecting Fig. || these discrepancies are not remedied by 
assuming a larger coupling ratio x. Furthermore, the sen- 
sitivity on the excited vector R\0) is not very large so that 
an explanation of the deviations in terms of modifications 
of the observables R in (||) is unlikely. Hence extensions 
of the model must indeed be considered. 

Numerous results favor the inclusion of a 4-spin cyclic 
exchange term -ffcyc as next important term [^l|-p4|. 
Especially in the spin ladder material La6Ca8Cu2404i 
about 10% cyclic exchange reconcile results for the spin 
gap, t he dispersion and the weighted spectral densities 
||g,||jlll. The term H^yc reads 



Hr 



2J< 



eye / , 



[(Si^iSi^i-i_i)(S2,iS' 



2A+1) 



(4) 



(Sl,iS2^i)(Si^i+lS2^i+l) — (Si_,;S2,i+l)(Si^i-)_iS2^i)] 



To assess the effect of a cyclic exchange term on the Ra- 
man hne shape quahtatively we include the 1-triplet dis- 



persion for various values of Xc 



f^cyc/ f^pc 



the 



inset in Fig. 0. The gap a;(7r) decreases on increasing 



f'cyc 



For Xcyc = 0.1 we find A = 0.337Jl in very 
good agreement with the result from neutron scattering 
A = 0.343J_L ^^. More important for the Raman line 
shape is that the dip in the dispersion Loik) at fc = 
is reduced on increasing Xcyc- The dispersion is mono- 
tonic for Xcyc — 0.10. Additionally, preliminary results 
on the 2-triplet terms show that the attractive inter- 
action between triplets is diminished by the cyclic ex- 
change. Thus we expect the pronounced zero inside the 
continuum of the Raman line to disappear on inclusion 
of cyclic exchange leading to an asymmetric broad peak 
with a shoulder at the high energy side as measured, see 
Fig. ^ and Refs. [qJIJi^I- We estimate that the deviation 
(i) is remedied by a 4-spin cyclic exchange term with 
Xcyc ~ 0.1 in agreement with other analyses ||30| , |35| , [ll| . 

The deviation (ii) in the sharpness of the main peak 
may also be remedied by the inclusion of the 4-spin cyclic 
exchange term. Alternatively, one might invoke a dop- 
ing effect, namely holes in the ladder, since the lines are 
considerably sharper for Sri4Cu2404i jlOll , in particular 
in [cc] polarization. Other results |1 l| , p6[[ 7 however, show 
that there are no holes in the ladders shedding doubt on 
an explanation in terms of holes in the ladders. 

In conclusion, we calculated the Raman response in 
the antiferromagnetic 2-leg S = 1/2 Heisenberg ladder in 
terms of dressed rung-triplets as elementary excitations. 
Our results are based on a continuous unitary transfor- 
mation introducing the number of elementary excitations 
as good quantum number. We demonstrated that the 2- 
triplet contributions dominate by far. Unexpectedly we 
found a 2-peak structure composed of a broad main peak 
and a secondary peak at higher energies for J|| ^ 0.8J^. 
From the inconsistency with the experimental finding for 
La6Ca8Cu2404i we infer that the model must be ex- 
tended for which a 4-spin cyclic exchange term is the 
best candidate to date. To exclude possible effects due 
to doping we urge for further investigations in undoped 
materials such as SrCu203. 
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